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Abstract 

The equilibrium quantum statistics of various anharmonic oscillator systems including 
relativistic systems is considered within the Wigner phase space formalism. For this purpose 
the Wigner series expansion for the partition function is generalised to include relativistic 
corrections. The new series for partition functions and all thermodynamic potentials yield 
quantum corrections in terms of powers of ft 2 and relativistic corrections given by Kelvin 
functions (modified Hankel funktions) K y (mc 1 2 /kT). As applications are treated the sym- 
metric Toda oscillator, isotonic and singular anharmonic oscillators and hindered rotators, 
i.e. oscillators with cosine potential. 

1 Introduction 

In recent years, the Wigner formalism as a phase space representation of quantum mechanics, 
quantum field theory and quantum statistics has found growing interest [1,2]. One of the main 
fields is the theory of anharmonic oscillator systems modelling various quantum systems, e g 
solitonic systems, quantum field theories, and transport processes in more complicated systems. 
The purpose of this paper is twofold, 

1. to generalize the Wigner series expansion of equilibrium phase space quantum statistics in 
order to include special relativistic systems, 

2. to show applicability and utility of the formalism by means of various examples. 

Up to now a small number of papers on relativistic quantum theory in the Wigner formalism 
have appeared for equilibrium as well as for nonequilibrium systems. The topics treated for 
nonequilibrium processes range from multiparticle production and kinetic theory [3] to cosmology 
[4]. In the realm of relativistic phase space quantum mechanics there exist only a few papers. 
Janussis et al. [5] starting from the Dirac Hamiltonian introduce a Wigner function with 4x4 
spinor components. Ali [6] discusses the quantization of relativistic systems on phase space whereas 
Kim and Wigner [7] describe a covariant phase space representation for harmonic oscillators. A 
relativistic Fermi gas is treated in the frame work of the grand canonical ensemble by Greiner et al. 
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[8]. A pseudo-classical phase space description of the relativistic electron in terms of Grassmann 
variables was given in several versions [9]. There are also papers on relativistic quantum kinematics 

l 10 ’ 11 ] 

Over the years, fundamental problems regarding the correct description of simple systems 
such as the classical relativistic harmonic oscillator and its quantization have been treated in a 
number of publications, see f.i. [12-15]. Important papers on the covariance problem are [7,16] ; 
a symplectic formulation of relativistic quantum mechanics has been given in [17,18]. We develop 
a Wigner formalism which is frame dependent. The focus of our attention is the comoving frame, 
i.e., we use the so-called synchronous gauge [19]. 

2 Relativistic Wigner Formalism 

One of the advantages of the phase space methods of quantum statistics - Wigner- Weyl formalism 
[l], coherent state methods [20], Husimi transforms [21] or Bargmann representation [22] - is the 
possibility to evaluate exactly the partition function by means of a phase space integration. In 
practice it is convenient to expand the phase space integral into the Wigner series in powers of h 2 ; 

[23] as a basis of semiclassical quantum^ statistic s. ;• . j 

The (non-relativistic) Wigner function can be defined as the Fourier transform of the j 


off-diagonal elements of the density operator (q | p | q 1 ) (for systems without spin) [l]: 

fw(q, P ) = ^ 7Ky J ^ ex P(*> z /*)<?-§ M?+ §> (0 

Here f is the number of the degrees of freedom of the system considered. The density operator is 
given in the coordinate representation by 

<« U I ?') = £ «9>n«9) (2) 

n 

with w n = exp (—PE n )/Z, where Hip„ = E n ip n , P — 1 fkT, and the partition function is 

Z = Tr exp (-0H) = £ exp(-0£ n ) (3) 

n 

Within the frame of the Wigner formalism the partition function can be expressed as 

Z-Jjdqdp^pj) , (4) 

where Clw = (exp(— 0H))w is the Wigner equivalent of the operator H = exp (—/3H) defined as 
SV(g,f>;/d) = J dzexp(xpz/h)( q - Z -\Cl\q+i) (5) 

The Bloch equation for the calculation of Clw within the Wigner formalisms reads 

= -Hw(q, P )co8(^A)Q w (q,pJ) , (6) 
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here Hw{q t p) is the Wigner equivalent of the Hamiltonian H t and A denotes the Poisson Bracket 
operator or symplectic differential operator 

a = ( 7 ) 

dp dq dq dp 

acting in (6) to the left and right hand sides. Up to now the formalism is completely exact. 
The partition function (4) can be calculated by means of a phase space integration without any 
knowledge of the energy levels of the system considered and without the necessity to sum up an 
infinite series of Boltzmann terms exp , the zero-point energy beeing already included. 

In order to elucidate the principle of our formalism, we consider a system with a single degree 
of freedom characterized by a Hamiltonian H = Ho{p) + U(q) with 

Ho(p) = 'J m 7 c* + p 2 c 2 — me 2 (8) 


and the potential energy V(q) . One has H w (q,P ) = H 0 (p) + V{q) . The Wigner series [17] for 
the partition function Z is given by 


Z 




$n(q,p;0) 


( 9 ) 


with $0 = 1,$*: n-th quantum correction. Solving the Bloch equation (6) for our relativistic 
quantum system yields 

*1 (q,pj) = V” + V* + H? V ") (10) 

where the primes denote differentiation with respect to the corresponding variables. We restrict 
ourselves to the first quantum correction. 

We will evaluate (9) with (10), this approximative partition function is denoted by Zoi ■ We 
begin with the p-integration: 

Zoi (tm) = Zo(tin) + Zl(*m) (H) 

with 2 0(t , n) = — J e' 6Bnlp) dp 

and Z 1( „ n) = V” + H? V"]dp 

Substituting pfmc = sinhu and using [43] 

gt coih u cos j, — K v {b) (12) 


where the K v denote the Kelvin functions (modified Hankel functions) one gets 

Z 0 (km) = 2^e +/W mcA' 1 (/9mc 3 ) 


(13) 
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and after some algebra 


» t^A,OW)-^Ay/W)+^jr s (/W)] + 

+ 2— e +flmc7 mch 2 — 0*V"<? x 
2irh 24 

X [-^Wmc’) + - -^K^mc^)] 


where the series expansions 

[i + _L_]- 3 / 2 = [1 + sinh 2 u]' 3/2 ~ 1 - | sinh 2 « + ^ sinh 4 « 
m 2 c 2 2 o 


( 14 ) 


(15) 


2 2 

— — [i + -£— l- 1 = sinh 2 u[l + sinh 2 u] _l ~ sinh 2 u - sinh 4 u (16) 

m 2 c 2 m 2 c 2 

have been used. 

Now the integration over the momentum space is done, and (13), (14) represent, as a main 
result, the partition function up to the first quantum correction, proportional to h 2 , and relativistic 
corrections up to second and fourth order in p/mc . To proceed further one has to take into 
account the potential V(q) of the system considered. Then the full partition function (in our 
approximation) is given by 


■'Ol 


- Z 0 + Zi = J dq exp(- / SV r (^))Zoi(*,„)(?) 


(17) 


Here we only mention that with 

*•<*> - ]fE e ~' 

in Zo(km) aIl< ^ 

- \fr, e ~‘ (1 + h ] • KM = \fE ‘" {1 + 3 i ) 

in Z X (km) the non-relativistic limit is reproduced correctly. 


K,{z) 



, 99 „ 

< 1 + K> 


3 Toda oscillator 

As a first example we choose the symmetric relativistic quantum Toda oscillator with the potential 

V(q) = V 0 (coshag — 1) (18) 

Its eigenvalues are not exactly known. If we compare (18) with the potential V(q) = mu>lq 2 / 2 

of the harmonic oscillator we find 

a = yj mull V 0 ( 19 ) 


326 


where w 0 is the (circular) frequency for harmonic (i.e. very small) oscillations. Toda molecules are 
discussed in [24], generalizations thereof in [25]. There exist a few papers on quantum statistics of 
toda chains [26] and Toda fields [27], The non-relativistic phase space quantum statistics of the 
symmetric Toda oscillator with (18) is treated in [28]. The dynamics of non quantum relativistic 
Toda lattices is the topic of [29-31]. 

The part Z 0 of (17) becomes 

Z 0 = 2 -i- me K^Pmc 2 ) e pv ° - K 0 { pV 0 ) (20) 

cl 

and the part Z\ can be expressed as 

Zi = 2 h 2 -ire e pmci Ki(0V o ) p 1 K 0 e™ x 

2 Th 

* _ + S Ars( ' 9racJ)1 + 

+ ^1- H (21 ) 

The formulae (20) and (21) yield in the nonrelativistic limit kT < me 2 the correct partition 
function of [28]. The evaluation of higher order relativistic and quantum corrections is straight 
forward. From the Wigner series for the partition function, corresponding series expressions for 
the thermodynamic potentials follow in the well-known manner. These results are published in 

[32]. 

4 Isotonic or Singular Oscillator 

Our next example is the quantum oscillator with the potential energy 

V(q) = V 0 [i - -) 2 (22) 

a q 

(normalized so that the potential minimum at g mm = a gives V(q mtn ) = 0). This asymmetric 
anharmonic oscillator appears in the literature under various headings in two different interpreta- 
tions: 

1. isotonic oscillator [33-36], radial oscillator (with centripetal barrier) [37-39] or inverse qua- 
dratic oscillator [40], 

2. singular oscillator [41] or nonpolynomial oscillator [42]. 

Interesting subject for oscillators of this type are the energy spectrum [38], canonical transfor- 
mations [37], ladder operators and coherent states [34-36], phase space path integrals and the 
dynamics and symplectic groups [37,41]. The connection with three-dimensional potentials was 
discussed thoroughly [39]. Here we add the second quantum correction to the partition function: 

(23) 
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where 


Q(q) = -VW + fi(2V'VW + ^-V" 2 ) - ^pV^Vll + ^70 3 V'* (24) 

£ O 

here is V' = dV(q)/dq , V^ 3 ) = V m , and so forth. For the isotonic oscillator potential (22) all 
integrals to be calculated are of the type (43): 


00 

/ 


,»-i e -(V+»r*) dx = K n/2 (2s/AB) 


(25) 


where the K -functions are again the modified Hankel functions (Kelvin functions) but here only 
of half integer index, n is an odd integer. These functions have a simple analytical shape, in 
particular one has 


K±\li,{ z ) = \J^ e 

K ±* d z ) = ][^ e ~'( 1 + b 

K ±sd z ) = ]ft e ~’( 1 + 7 + 7> ) 


(26) 


In (25) one has for the isotonic oscillator 

A = fiVo/a 3 , B = pV 0 a 3 } z = 2>/AB = 2pV 0 . 
The exact classical (i.e. semiclassical ) partition function becomes 

Z cl = kT/hui 


where 


w = >/8Vo/ma J 


(27) 


(28) 


(29) 


is the classical (angular) frequency of our oscillator in the limit of very small (i.e. harmonic) 
oscillations. The partition function (28) is formally the same as the well-known expression for the 
classical harmonic oscillator with (angular) frequency w. 

A simple but somewhat lengthy calculation then yields the first two terms of the partition 
function 

kT . 1,JW, 


M r t x / v 0/ „ 3 kT 


(30) 


where the expression in the square brackets being the famous h 2 / 24 quantum correction. The 
higher order quantum corrections can be calculated in a similar way using (25). Again all integrals 
lead only to the handy Kelvin functions with half-integer index. As is well known one gets 
the thermodynamic functions for a system of N oscillators from the partition function. The 


328 


(Helmholtz) free energy (here identical to the free enthalpy G) , the internal energy U(here identical 
to the enthalpy H), the entropy S and the heat capacity C = dU/dT are given, respectively, by 


F = -NkT In Z 
U = NkT 2 ^ In Z 

S = Nk{\n Z + T ^ In Z) 

c = ^^( T3 ^ ,nZ ) • (31) 

Starting from (30) these functions can be written in the form of classical contribution plus quantum 
correction up to terms with h 2 using the series expansion of the logarithm function: 

' - <£>'< i+ !£» (32) 
* - (33) 

S = Nk In + Nk{l + ( 34 ) 

c = < 35 ) 


More detailed estimates show that these semiclassical expressions are of considerable accuracy for 
atomic oscillators above T = 100 K, the errors being only a few percent and enlarging in the order 

F,U,S,C. 

5 Hindered Rotators 

The internal rotation of molecular groups (f.i. the methyl group) within molecules (f.i. ethane, 
polymers) or molecular complexes (f.i. C^Hs OHi ) is a topic of large interest since many years 
[44,45]. Related motions are the rotation of molecules adsorbed on solid surfaces [46,47] and 
the torsional oscillations of molecules or ions in the interior of molecular crystals [48,49] The 
contributions of the corresponding degrees of freedom to the thermodynamic functions and to the 
equilibrium and rate constants are, generally, important. In calculations in the frame of statistical 
thermodynamics these rotators were mostly treated as (classical) free rotators [50] or sometimes 
as (harmonic) torsional oscillators [51]. 

However the rotators have to be classified, in most cases, as hindered rotators or strongly 
anharmonic torsional oscillators, the more the lower the temperature is. There exists only a few 
number of papers on the quantum statistics of hindered rotators [52-56]whereas the literature 
on the quantum statistics of various types of anharmonic oscillators is abundant. The reason is 
that the Schrodinger equation for a hindered rotator is of the Mathieu-Hill type whose energy 
eigenvalues can be calculated only approximately. Therefore, the evaluation of the partition 
function and the thermodynamic functions is possible only with numerical methods; closed general 
formulae cannot be derived from this approach [52,53]. Here , the Wigner formalism shows its 
advantages. 
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The main quantities characterising a hindered rotator are the shape and the strength of the 
potential barrier. For a symmetric rotator, the hindering potential can be given by 

V(4>) = ^Vo (1 - cosndO (36) 

where n is the symmetry number (f.i. n = 3 for the methyl group) or, more generally, by Fourier 
series [57] 

V(<f>) = ^(1 - ^ a k coskn<f>) (37) 

* 

In the most simple case V{4>) is dominated by the repulsive interaction with a single near-neighbour 
atom within the molecule (or crystal) considered; then all coefficients a* for k > 2 are zero. The 
presumption is also that the cosine series converges rather rapidly. A reasonable representation of 
the barrier can then made with the first two terms. For the important case of n = 3 one has 

V(4>) = -[V 3 (l — cos3d>) — V® cos 6 ^] (38) 

2 

The effect of adding a small sixfold potential Vgcos 6<f> to the threefold main potential V 3 cos 3 <£ 
is to broaden the potential wells at the expanse of the barriers. (If V$ is negative then the wells 
are narrowed and the barriers broadened.) For many acetates and other molecules and complexes 
containing the methyl group we have < V 3 (f.i. Vjs/Va as 0.15 for acetates) [57]. 

For a hindered rotator with a simple cosine potential the canonically conjugate coordinates 
are given by the rotation angle <f> and the angular momentum L } respectively; and the mass is to 
be replaced by the moment of inertia I. 

In order to have a convenient comparison with the limiting case of a harmonic oscillator 
(Vq — > 00 ) we introduce the angular frequency Cl of of the harmonically oscillating torsional 
vibrator by 



With (36) the classical partition function is written a s 

Zd = j e ICO,n *d<i> (40) 

0 

where x = /?V"o/2. We get the well-known expression 


Z d = Z,e- z I 0 {x) (41) 

where 

n n 

is the classical is the classical partition function of a free internal rotator of symmetry number n, 
and Io(x) is the modified Bessel function (or the Bessel function Jo(ix) of the purely imaginary 
argument ix) of order zero [43]. 
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(43) 


The evaluation of the quantum corrections Z(h 3 ) and Z(h 4 ) yields with u — at first 


Z(h 3 ) = Z } e-*^ J e rco, "*(-^cosnd> + sin 2 n<f>) dcj> 
0 

4 2T 

Z(h*) = Z f e~ z ^ J e IC0,n * ^(icosn</> — 2sin 2 n<^ + 

0 

3 11 5 

+ - cos 2 n<f> -x sin 2 nif> cos n<f> + — x 3 sin 4 n<f > ) d<f> 

2 6 24 


(44) 


All integrals in (43) and (44) can be expressed by Bessel functions J v {z ) of purely imaginary 
argument z = »£, where v is a positive integer index [43], 


2t 


. h 


o 

or by the modified Bessel functions 


e »xcoi n<f> 


h(x) = \~ v J v (ix) 


One gets 


1 


(45) 


(46) 


(47) 


Z(h 3 ) = -Zfe-'-pXhufhix) 

Z(h 4 ) = +Zfe- I ^0Xhw) i {-jh(x) + l[I o (x)-jh(x)) - 

- ^[A(*) -/,(*)] + |/ 3 (*)} . 

By means of the recursive relations [43] 

h-i(x) - I*+i(x) = 2ul v (x)fx 

we have finally the total partition function as a compact expression 

- - -,mn ewhw nw 2A(x) 

Z-Z,e /«(i){1 24 + 5760 ( ( /„(*)" 

The factor in curly brackets represents the quantum correction to the classical partition function 

(4i). 

In the case of a modulated cosine potential given by 


(48) 


(49) 


(50) 


V(4>) = -Vb(l - cosn</>- ccos2 ntf) 

2 


(51) 
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where e » 0.1 for most molecules or molecuar complexes, at first we calculate the classical partition 
function written in the shape 


Zd 



oo p p 

e r co. n# + £ - cos*' 2n<j>) d<t> 

V-\ 


(52) 


With [43] 

fj — 1 y v 

cos lM_1 2 n<f> = Y^, [ k 1 J cos (2^i — 2k — l)2n<£ (53) 

Jfc= 0 ' ' 

cos J m 2 n<f> = 2 cos2(#i- k)2n<f> + Q^)} (54) 

one gets for Z c i a series expansion in terms of Bessel functions /j( a„- 3 *-i)(i) and A( m _*)(i), 
respectively. Clearly this series converges rather rapidly, and in practice one needs only terms up 
to the order e or e 7 . Therefore we use the approximation 


2t 
1 t 

Z c i = Z/e _I — / e IC0,n * [i + C x cos2n0 + e 3 — cos 3 2n$ d<f> 
2 x J 2 


(55) 


and get, by use of (45), after elementary transformations 


(56) 


Comparing this expression with the classical partition function (41) without the potential propor- 
tional to ecos2n<£> we identify the factor in square brackets in (56) as the potential modulation 
correction factor in Z c i . With the help of the recursive relations (49) it is again possible to express 
the whole correction factor in terms of the Bessel functions /o(0 and I\(() only. 

We content ourselves with the first quantum correction Z(h 2 ) and write 


Z(h 7 ) = Z f e- z ± j e *<" n +[l+'R l + e 7 R i ]^{-V"+^V' 7 }d<t> 
0 

where 

R\ = x cos 2 n<f> 

x 2 x 2 

R 2 = — cos 2 2 n<p = —(1 + cos 4n<£) 

i 

V' = -V 0 (nsinnd> + e2nsm2n<f>) 
z 

V" = -Vo(n 3 cosn$ + f4n 3 cos2n^) 

2 


(57) 


(58) 

(59) 

(60) 
(61) 
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After somewhat lengthy calculations using elementary formulae for trigonometric functions of 
multiple arcs and (45), (46), one ends with 


with 


Z(h') = + «(§/.(*) + «*)) - c'«} 


a = x(6Iq(x) - -/i(s) + 4/ 2 (x) - / 3 (x) - 2/ 4 (x)) + x 7 (Iz(x) - /s(x)) 


(62) 


(63) 


This quantum correction again contains terms proportional to e and c 3 due to the potential 
modulation. From now we neglect the term e 2 a. For the free energy of this simple hindered 
rotator we get: 

F = F cl + F{h 3 ) (64) 


where 

F cl = -FT In Z f - RT { In J 0 (x) - x) 
is the free energy of the semiclassical approximation [56], and 

r ,^ (^u/^/xCx) 

F(h)-KT ■ 


(65) 


( 66 ) 


Now we come to the hindered rotator with the modulated cosine potential. First we write the 
full partition function according to (56) 

Z = Z cl + Z{h 2 ) 


„ .... , .(teW'/tt!) ,3/,(l) /,(*) 

-z ;e + z ^)] 




(67) 


where A,B are temporary abbreviations ( notice that Jo(x) > 1 , x = £Vb/2). The free energy 
now becomes 


F = -ET ln(Z/e-*/ 0 ) - RT{A - - \{A 2 - ^~^AB)} 


( 68 ) 


Here again we take the quantum correction up to terms in h 3 and consider the potential modulation 
correction up to terms in e 2 in the semiclassical contributions and up to terms in e within the 
quantum correction. With these approximations we get finally: 


H*) 




F = -Iff In Z; -MXIn W-. )-J*n«x^ + ^l + ^ -3(^)1 + 

, pr P&fl r hW I gT r 3 Mfl ■ Mil 

24 l / 0 (x) + n 2/o(x) + / 0 (x) I2(x) }l 


(69) 
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